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Abstract. An ideal quantum walk transitions from one vertex to another with perfect 
fidelity, but in physical systems, the particle may be hindered by potential energy 
barriers. Then the particle has some amplitude of tunneling through the barriers, 
and some amplitude of staying put. We investigate the algorithmic consequence of 
such barriers for the quantum walk formulation of Grover’s algorithm. We prove 
that the failure amplitude must scale as 0{l/\fN) for search to retain its quantum 
0{'/N) runtime; otherwise, it searches in classical 0{N) time. Thus searching larger 
“databases” requires increasingly reliable hop operations or error correction. This 
condition holds for both discrete- and continuous-time quantum walks. 


PACS numbers: 03.67.Ac, 03.67.Lx 

1. Introduction 

Classical random walks, or Markov chains, serve as the foundation of many classical 
probabilistic algorithms [1]. It is no surprise, then, that their quantum analogues, 
quantum walks similarly form the basis of many important quantum algorithms. 

This includes quantum search [5], element distinctness [6], and evaluating NAND trees 
[7], all of which provably yield polynomial speedups over the best possible classical 
algorithms. In other applications Pi, the speedup can even be exponential, and any 
quantum algorithm can be efficiently simulated by a quantum walk [lOj . 

These algorithms assume that the quantum particle transitions from one vertex 
of a graph to another without hindrance. In physical systems, however, this may not 
be the case. As explained in various physical processes can prevent the particle 
from hopping, and the obstacles can be modeled as potential energy barriers. Then 
the randomly walking quantum particle must tunnel through the barriers in order to 
transition between vertices. As shown in du, even with the barriers, a quantum walk 
on the one-dimensional line still exhibits the ballistic propagation that is characteristic 
of quantum walks [1^. Despite this, the work leaves open whether the barriers have 
algorithmic consequences. 

In this paper, we investigate whether potential barriers have algorithmic 
consequences by exploring their effect on the discrete- and continuous-time quantum 
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walk formulations of Grover’s algorithm na. In the next section, we introduce 
discrete-time quantum walks and apply them to the unstructured search problem with 
potential barriers. We prove that, for the search to keep its quantum 0{'/N) runtime, 
the amplitude of failing to hop must go to zero sufficiently quickly, specifically as 
0{1/\/N). Otherwise, the algorithm searches in 0{N) time, losing its speedup over 
classical. Thus the potential barriers have significant algorithmic consequences; either 
the potential barriers must be asymptotically eliminated, or error correction [T3] must 
be utilized. Following, we analyze search with continuous-time quantum walks with 
potential barriers, arriving at the same condition for the error amplitude. Finally, we 
end with concluding remarks. 

2. Discrete-Time Quantum Walks 

Grover’s algorithm solves the unstructured search problem, which is equivalent to 
searching the complete graph of N vertices for a marked vertex [ladS], as illustrated 
in figure [ 1 } First let us introduce the quantum walk without the searching component. 
The vertices of the graph label computational basis states {|1), |2),..., |A^)} of an N- 
dimensional “vertex” Hilbert space C^. In discrete-time, these vertices only support a 
trivial quantum walk [HI EH], so assuming the graph is d-regular, we define an additional 
d-dimensional “coin” Hilbert space C'’* spanned by the d directions in which the particle 
can hop from one vertex to another. So the system evolves in For the complete 

graph, each vertex is connected to the N — 1 other vertices, so d = — 1. 

Without potential barriers, each step of the quantum walk is performed by applying 
a “coin flip” followed by a hop/shift: 

Uo = S-{In®Co), 

where Cq is the “Grover diffusion” coin [5], 

Cq = 2|Sc)(Sc| ~ hdy 

where |sc) = N)/v^ is the equal superposition over the coin space, and S is the 

flip-flop shift [15] that causes the particle to jump from one vertex to another and then 
turn around {e.g., S']!) G |1 —)■ 2) = |2) ® |2 —)■ 1)). 

Now with potential barriers, the particle has some amplitude a of successfully 
hopping and some (related) amplitude (3 of staying put, leading to the following faulty 
shift; 

S ^aS + pi. 

Since S is both unitary and Hermitian, aS + [31 is unitary provided lap -|- |/lp = 1 and 
a[3* + (3a* = 0. Given these constraints, we parameterize 

a = cos 0 and /I = i sin <p. 

Then the quantum walk operator is 

Uq = (cos(0) S + isin(0) J) ■ (/at ® Cq). 
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Figure 1. The complete graph with N = 6 vertices. A vertex is marked, as indicated 
by a double circle. Identically evolving vertices are identically colored and labeled. 


Note that 0 can be interpreted as a failnre parameter. When 0 = 0, the walk is ideal 
and there are no potential barriers. Bnt as 0 —)■ 7r/2, the barriers get worse and worse, 
cansing the hop to fail more. 

If the particle begins in the nniform distribntion, i.e., |0o) = k,;) ® |sc), where 
Is.;;) = \i)/^/N is the eqnal snperposition over the vertex space, then the potential 

barriers makes no difference, i.e., f/o|0o) = |0o) and t/o|0o) = |0o)- To tnrn this 
qnantnm walk into a search problem, we look for a particnlar “marked” vertex |a), as 
shown in hgnre by qnerying an oracle Ra that flips the phase of the marked vertex, 
i.e., Ra\a) = —|a) and Ra\x) = |a;),Va; ^ a. Qnerying the oracle with each step of the 
qnantnm walk [19], the search is performed by repeatedly applying 

U = (cos(0) S + isin(0) I) • {1^ 0 Cq) • {Ra 0 h) (1) 

to the initial eqnal snperposition |0o) = l-s,,) 0 |sc). As proved in [20], when the particle 
freely transitions from one vertex to another withont the potential barriers {i.e., 0 = 0), 
the snccess probability reaches 1/2 after 'k\/N/2'/2 applications of U, as shown in 
hgnre 1^ which resnlts in a Q{y/N) search algorithm with the expected constant nnmber 
of classical repetitions to boost the snccess probability near 1. As in typical discrete¬ 
time qnantnm walk algorithms, this snccess probability of 1/2 assnmes that only the 
position of the particle is measnred. If the internal state of the particle is additionally 
measnred, however, the snccess probability can be boosted to 1 IZB. 

With the potential barriers {i.e., 0 > 0), the particle has some amplitnde of not 
hopping. Even for small valnes of 0, this can signihcantly impair the search, also shown 
in hgnre Let ns hnd how large 0 can be snch that we still search in 0 (a/]V) time. In 
doing so, we will get a sense for how small the potential barriers mnst be for qnantnm 
walks to search qnickly. 

To do this, we explicitly work ont the evolntion of the system. As shown in hgnre 
there are only two types of vertices: the marked red a vertex and the nnmarked white 
b vertices. By symmetry, the b vertices evolve identically. Now the a vertex can only 
point towards b vertices, while b vertices can either point towards the a vertex or other 
b vertices. Thns the system evolves in a 3D snbspace spanned by these vertices and 
directions: 


1 
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Figure 2. Success probability as a function of the number of applications of U 
for search on the complete graph with N = 1024 vertices with (j) = 0 (no potential 
barriers), = 0.02, and (j) = 0-04 corresponding to the solid black, dashed red, and 
dotted green curves, respectively. 
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The (unnormalized) eigenvectors and corresponding eigenvalues of this are 
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/ — cot I \ 


cot I 


V 




/ 


and 
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( (1-e-‘(^'"+'^))csc0 \ 

cos 9 + CSC 9 
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,±icr 


where 


cos a = 


(1 + cos 6*) COS0 


sin a = 


+ cos 9Y cos^ (j) 


2 ’ 2 
It will be useful in our analysis to have the sum and difference of "0+ and '0_ 
because they evolve to each other up to an overall factor of i (he., up to a phase) in 
7r/2cT applications of the search operator U ([^. That is, 

± ^l)_) = (e^'")^/2<7^+ ± = i{'ip+ T Y-)- (2) 

In particular, the sum and difference of 'ip+ and Y- are 

^ 2(1 — e“''^cos(T) csc^ \ 
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Pluging in for a and 9, they become 

Y+ + Y- = 
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We can use these sum and difference formulas to find the evolution of the search 
algorithm when 0 scales less than, equal to, or greater than l/\fN. To simplify the 
calculation, we substitue (f) = cj\fN and equivalently consider when c scales less than, 
equal to, or greater than a constant. Then Taylor expanding for large N, we get 

i ‘"+o(w) ) 
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V’+ - '0- = 




(4) 


V 0 / 

Now let us consider each scaling of 0 = c/ \/iV or c as separate cases. 

Case 1: (j) = o{ 1/\/N) or c = o(l). When 0 scales less than l/\fN, or equivalently 
when c scales less than a constant, the sum ([^ and difference Q of and xjj- have 
leading-order terms 

/o\ ( 


V'+ + V'- 
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VV 





Then the initial state is approximately 
iV'o) ~ \bh) ^ ^ + . 

Using after 7r/2cr applications of [/, the system evolves to 
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which is half in \ab) and half in \ba). This corresponds to a runtime of 
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So if 0 scales less than l/y/N, the system evolves from I^q) ~ 1^^) to being half in \ab) 
and half in \ba) after tt\/N /2\/2 applications of U , which results in a success probability 
of 1/2 when measuring the position of the particle due to the \ab) piece. This result is the 
same as without the potential barriers [20], and so the potential barriers are too small 
to affect the search for large N. An example of this is shown in hgurej^with 0 = 
and N = 1024, 4096, and 16384. For each of these, the success probability reaches 1/2 
at respective times 7r\/l024/2\/2 36, 7r\/4096/2-\/2 th 71, and 7rA/l6384/2-\/2 th 142, 

as expected. 

Case 2: 0 = 0(l/\//V) or c = 0(1). When 0 = c/\/N with coefficient c constant, 
the sum ([^ and difference (|^ of 0+ and 0_ have leading-order terms 
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Figure 3. Success probability as a function of the number of applications of U for 
search on the complete graph with (j) = 1/N^P and N = 1024, 4096, and 16384 vertices 
corresponding to the solid black, dashed red, and dotted green curves, respectively. 


Since the initial state |'0o) ~ 1^^) is ronghly half this, after ti/ 2a applications of U, the 
state of the system is (§ 
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“ 2 
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Measnring the position of the particle yields a snccess probability given by the sqnare 
of the hrst term, which corresponds to \ab): 

-2i 
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The corresponding rnntime is 

TT TT TT TT 

^ _ 

2a 2sin(T ^/A(|P + (3 + cos6l)(l — cos 6*) + 8/N 

TT 

?T2 

An example of this is shown in hgnre with c = 1 and N = 1024, 4096, and 16384. 
Each of these reach a snccess probability of 1/(1^ + 2) = 1/3 at respective times 
7rVl024/2Vl2 + 2 29, Ti^/Wm/2^/l‘^ + 2 58, and 7rVl6384/2Vl2 + 2 ^ 116, as 

expected. While the nnmber of steps and snccess probability are smaller than in the 
hrst case, their scalings are nnchanged, so we still achieve Grover’s 0(a//V) steps with 
snch potential barriers. 

Case 3: 0 = ijj(l/\fN) or c = a;(l) (and </> = o(l) or c = o{\/N)). When cj) scales 
larger than l/y/N (bnt still less than a constant), or eqnivalently when c scales greater 
than a constant (bnt less than y/N), the snm and difference of "0+ and '0_ for large N 
are 


V’+ + V’- = 

/ iy/N(j) \ 

+ 
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/-iy/N(j)\ 


[ 2 J 
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Figure 4. Success probability as a function of the number of applications of U for 
search on the complete graph with (ft = l/\fN and N = 1024, 4096, and 16384 vertices 
corresponding to the solid black, dashed red, and dotted green curves, respectively. 


Then 


2V’+ = (^+ + V'-) + (V'+ - -0-) ~ I 0 


= 2\hh) 

Thus the system approximately begins in an eigenstate: 


0^ 

/ 


iV’o) ~ \hh) ^ ^/^+, 

so it fails to evolve (apart from a global, unobservable phase), and hence the search 
fails. That is, measuring the state at later time equates to measuring the initial equal 
superposition state, which gives the marked vertex with probability 1/iV, which is 
equivalent to classically guessing and checking. An example of this failure is shown 
in hgure with (j) = 1/N^^'^. As N increases, the success probability evolves less and 
less from its initial value of 1/N. 

Case 4- <P = 0(1) or c = Q{\/N). When 0 scales as a constant, which is the largest 
it can scale since 0 = 7r/2 corresponds to the potential barriers stopping all transitions, 
then the behavior from Case 3 persists—for large N, the system approximately begins in 
an eigenstate, and so the success probability does not evolve. This is shown in hgure 
where for constant 0, increasing N causes the success probability to evolve less and less. 

Thus we see an abrupt change in the behavior of the algorithm for large N, 
depending on 0: when 0 = 0(l/\/iV), the full quantum quadratic speedup is achieved, 
and when 0 = no speedup over classical is provided. This is summarized in 

table Since the faulty shift is aS + f3I with a = cos0 and 0 = isin0, the amplitude of 
the particle failing to hop is 0 = i sin 0 i0 for small 0. Therefore the failure amplitude 
0 must scale as 0(1/a//V) for search to retain its quantum speedup; otherwise, the 
algorithm searches in classical time. This means the failure amplitude must decrease as 
N increases, so searching larger “databases” requires increasingly reliable hop operators 
in the absense of error correction. In follow-up work with Ambainis [H], we show how 
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Figure 5. Success probability as a function of the number of applications of U for 
search on the complete graph with 4> = 1/N^P and N = 1024, 4096, and 16384 vertices 
corresponding to the solid black, dashed red, and dotted green curves, respectively. 



Figure 6. Success probability as a function of the number of applications of U for 
search on the complete graph with </> = 0.02 and N = 1024, 4096, 16384, and 65536 
vertices corresponding to the solid black, dashed red, dotted green, and dot-dashed 
blue curves, respectively. 


the coin and oracle operators can be modified to partially offset the potential barriers, 
recovering the 0{\fN) runtime so long as 0 does not approach 7r/2. To do so, we 
novelly interpret the quantum walk algorithm as an amplitude amplification algorithm, 
then adjust the phases applied to boost the success probability [221 1^ - 

3. Continuous-Time Quantum Walks 

We end by discussing continuous-time quantum walks, which do not require the 
additional “coin” space, so the system evolves in the vertex Hilbert space C^. The 
system begins in the equal superposition over the vertices |s^), and without potential 
barriers, evolves by Schrodinger’s equation with Hamiltonian 

H = —yH — |a)(a|. 
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Table 1. Summary of discrete-time quantum walk search on the complete graph with 
potential barriers for large N. In the last case, the system stays in its initial uniform 
state, so the quantum algorithm is equivalent to classically guessing and checking. 


Barrier 

Runtime 

Success Probability 

Example (s) 

d = oAls/N) 


1 

2 

Figure 

3 


ij) = c/y/N, constant c 


1 

c2+2 

Figure 

1 


(j) = loA/'/N) 

Not Applicable 

1 

N 

Figures [^c 

in 

dg 


where 7 is an adjustable parameter corresponding to the jumping rate (amplitude per 
time), and A is the adjacency matrix of the graph = 1 if vertices i and j are 
adjacent). For the complete graph, all the vertices are connected to each other, so the 
adjacency matrix is 


/O 1 ... 1\ 
1 0 ... 1 

Vi 1 ... oy 


and it effects the quantum walk ng. Note since the complete graph is regular, using the 
adjacency matrix is equivalent to using the graph Laplacian m- As shown in figure 
there are only two types of vertices a and b, so we can write H in the 2D subspace 
spanned by |a) and \b) = \x)/\/N — 1: 


FT = -7 




x/aFi\ 
N-2 j 


As shown in [25], this has eigenstates 


IFi) 0^ + 


1 — 'jN ± AE 
27 vF 


a) 


with gap in the corresponding eigenvalues Eq and Ei 


AE = El-E q = V(1 - 7 FF 47 . 


When '^N takes its critical value of 1, the energy gap is AE = 2/vF with eigenstates 
1 - 00 , 1 ) cc |s^) ± |a), so the system evolves from |st,) to |a) with probability 1 in time 
tt/AE = 7rvF/2 [ini 123 120]. This behavior is retained near the critical 7 when 
'-)N = 1 -|- o(1/a//V) for large N. When '-^N = 1 -|- c/N for constant c, we get 
AE = ^J{(A + A)JN and I'lpo,!) oc |s^) + (—c ± \J(A E 4)/2|a) = Is^,) + o(l)|a) for 
large Af, so the system evolves from |s^) to |a) with constant probability in time 
T\\fN (A P 4. Finally, when y/F = 1 -f a;(l/vF), we get that |' 0 o) or iV'i) equals 
|s.u) for large Af, depending on if the deviation is positive or negative, respectively. So 
the system begins in an eigenstate and fails to evolve beyond acquiring an unobservable 
phase, which is consistent with degenerate perturbation theory [20]. Thus we require 
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Figure 7. Success probability as a function of time for continuous-time search on the 
complete graph with N = 1024 vertices, 7 = 1/A^, and e = 0 (no potential barriers), 
e = 0.02, and e = 0.04 corresponding to the solid black, dashed red, and dotted green 
curves, respectively. 


'yN = 1 -b 0{1/\/N) to achieve a 0(\/iV) search, and otherwise the search is no better 
than classically gnessing and checking. 

Now with potential barriers, the amplitnde that the particle hops from one vertex 
to another is decreased, say by e, and the amplitnde of doing nothing is increased by 
the same amonnt. This effectively modihes the adjacency matrix to be 


1-e 

1 - e (iV- l)e 


1-e \ 

1 -e 


V 1 -e 1 -e ... 

Bnt this is simply [N — l)eJ + (1 — e)A, and since the hrst term is a mnltiple of the 
identity matrix, it yields no observable effect and can be dropped. Thns the search 
Hamiltonian is effectively 


H = — 7(1 — e)A — |a)(a|. 

The actnal critical 7 for this walk solves 7 iV(l — e) = 1, which implies 'yN = 1/(1 — e) ~ 
1 -b e for e = o(l). If 'yN is still chosen to be its barrier-free valne of 1, then we are e 
away from the trne valne. Bnt above, we showed search in 0(a/]V) time is still possible 
when e = (9(l/\/iV), which is the same constraint as the discrete-time qnantnm walk. 
That the system evolves less and less as the error e grows is illustrated in hgure 


4. Conclusion 

We have shown the effect of potential barriers hindering a randomly walking quantum 
particle from searching on the complete graph of N vertices. In discrete-time, the 
amplitude of not hopping must scale less than or equal to l/y/N for the search to 
achieve Grover’s Q{\/N) runtime without error correction. Otherwise, no improvement 
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over classical is achieved. Similar behavior holds for search by continuous-time 
quantum walk. So even though a quantum walk on the one-dimensional line retains its 
characteristic ballistic dispersion with potential barriers HU, our result indicates that 
the barriers can have signihcant algorithmic consequences. Further research includes 
how the barriers affect other algorithms based on quantum walks, and the effects of 
non-uniform or random barriers. 

Acknowledgments 

This work was supported by the European Union Seventh Framework Programme 
(FP7/2007-2013) under the QALGO (Grant Agreement No. 600700) project, and the 
ERG Advanced Grant MQG. 

References 

[1] Norris J R 1998 Markov Chains Cambridge Series in Statistical and Probabilistic Mathematics 

(Cambridge University Press) ISBN 9780521633963 

[2] Aharonov Y, Davidovich L and Zagury N 1993 Quantum random walks Phys. Rev. A 48 1687-1690 

[3] Ambainis A 2003 Quantum walks and their algorithmic applications Int. J. Quantum Inf. 01(04) 

507-518 

[4] Kempe J 2003 Quantum random walks: An introductory overview Contemp. Phys. 44(4) 307-327 

[5] Shenvi N, Kempe J and Whaley K B 2003 Quantum random-walk search algorithm Phys. Rev. A 

67 052307 

[6] Ambainis A 2004 Quantum walk algorithm for element distinctness Proceedings of the 45th Annual 

IEEE Symposium on Foundations of Computer Science FOCS ’04 (IEEE Computer Society) pp. 
22-31 

[7] Farhi E, Goldstone J and Gutmann S 2008 A quantum algorithm for the Hamiltonian NAND tree 

Theory Comput. 4(8) 169-190 

[8] Childs A M, Cleve R, Deotto E, Farhi E, Gutmann S and Spielman D A 2003 Exponential 

algorithmic speedup by a quantum walk Proceedings of the 35th Annual ACM Symposium on 
Theory of Computing STOC ’03 (New York, NY, USA: AGM) pp. 59-68 

[9] Childs A M, Schulman L J and Vazirani U V 2007 Quantum algorithms for hidden nonlinear 

structures Proceedings of the 48th Annual IEEE Symposium on Foundations of Computer Science 
EOCS ’07 pp. 395-404 

[10] Childs A M, Gosset D and Webb Z 2013 Universal computation by multiparticle quantum walk 

Science 339(6121) 791-794 

[11] Wong T G 2015 Quantum walk on the line through potential barriers Quantum Inf. Process. 15(2) 

675-688 

[12] Ambainis A, Bach E, Nayak A, Vishwanath A and Watrous J 2001 One-dimensional quantum 

walks Proceedings of the 33rd Annual ACM Symposium on Theory of Computing STOC ’01 
(New York, NY, USA: ACM) pp. 37-49 

[13] Grover L K 1996 A fast quantum mechanical algorithm for database search Proceedings of the 28th 

Annual ACM Symposium on Theory of Computing STOG ’96 (New York, NY, USA: ACM) pp. 
212-219 

[14] Ambainis A and Wong T G 2015 Correcting for potential barriers in quantum walk search Quantum 

Inf Comput. 15(15&16) 1365-1372 

[15] Ambainis A, Kempe J and Rivosh A 2005 Coins make quantum walks faster Proceedings of the 

16th Annual ACM-SIAM Symposium on Discrete Algorithms SODA ’05 (Philadelphia, PA, USA: 
SIAM) pp. 1099-1108 


Quantum Walk Search through Potential Barriers 


13 


[16] Childs A M and Goldstone J 2004 Spatial search by quantum walk Phys. Rev. A 70 022314 

[17] Meyer D A 1996 From quantum cellular automata to quantum lattice gases J. Stat. Phys. 85(5-6) 

551-574 

[18] Meyer D A 1996 On the absence of homogeneous scalar unitary cellular automata Phys. Lett. A 

223(5) 337-340 

[19] Wong T G and Ambainis A 2015 Quantum search with multiple walk steps per oracle query Phys. 

Rev. A 92 022338 

[20] Wong T G 2015 Grover search with lackadaisical quantum walks J. Phys. A: Math. Theor. 48(43) 

435304 

[21] Prusis K, Vihrovs J and Wong T G 2015 Doubling the success of quantum walk search using 

internal-state measurements {a}rXiv: 1511.03865 [quant-ph] 

[22] Long G L, Li Y S, Zhang W L and Niu L 1999 Phase matching in quantum searching Phys. Lett. 

A 262(1) 27-34 

[23] Hpyer P 2000 Arbitrary phases in quantum amplitude amplification Phys. Rev. A 62 052304 

[24] Wong T G, Tarrataca L and Nahimov N 2016 Laplacian versus adjacency matrix in quantum walk 

search Quantum Inf. Process, pp. 1-20 {a}rXiv; 1512.05554[quaint-ph] 

[25] Wong T G 2014 Nonlinear Quantum Search PhD dissertation 

[26] Janmark J, Meyer D A and Wong T G 2014 Global symmetry is unnecessary for fast quantum 

search Phys. Rev. Lett. 112 210502 



